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1 . INTRODUCTION 


/ 3 » 


The problems Involved in studying the earth’s natural 
resources, just as any other problem, are connected with the 
processing of large files of information needed primarily in the 
visual representation of this information. The quasi-photo- 
graphy of video information obtained on the alphanumeric printer 
of a computer [1] gives the researcher a convenient opportunity 
not only to evaluate the content of the Incoming information, 
but also to identify the data recorded on magnetic tape and the 
image of the earth's surface with a precision down to the element 
of resolution. 

However, even at the very beginning stage of research, it 
is necessary to study the material numerically. Such a study is 
impossible without an extensive use of computer technology. The 
construction of graphs, histograms, digital print-outs, compu- 
tation of elementary statistical characteristics, etc. — this is 
far from a complete enumeration of the problems which the 
researcher has to know how to solve when processing large data 
files. 

The data recorded on magnetic tape from remote sensing of 
the terrestrial surface are color vectors assembled in a line [2] . 
The number of vectors in a line can vary — this depends on the 
properties of the recording apparatus. Each sequence of lines 
recorded on tape is an image. 

Each color vector Q ,j l (°\ • • > a *) can have up to eignt compon- 
ents which characterize the intensity of the electromagnetic 
radiation reflected from an element of resolution in one of the 
spectral ranges. During recording of the radiation, the signal A 
is modulated in such a way that each of these components a s can 
have a maximum value of 255 nominal units, i.e., it can be placed 
in a single byte. 

"Numbers in margin indicate foreign pagination. 
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The position of an element of resolution from an image is 
determined by the number of the line and the number of the color 
vector in the line. Thus the entire image can be regarded as a 
matrix whose elements are the color vectors °ij 

The proposed program makes it possible to obtain a chart of 
the change in the signal (along a line or a column of the image) 
in digital or graphic form; to calculate the vector magnitudes 
of the mathematical expectation and dispersion; to construct 
histograms of the distribution of the signal in each spectral 
range. The program consists of several blocks written in the 
languages ASSEMBLER and FORTRAN for the Disk Operating System 
of electronic computers of the Unified Series (AOC/EC) [3] . The 
reading of magnetic tape (ufiT , storage, rearrangement of files 
and the development of charts is carried out by blocks written 
in the ASSEMBLER language. This makes possible more complete 
use of the machine’s capacity. The storage files are transferred 
for subsequent processing by FORTRAN modules, which facilitates 
programing processing. 

Preliminary processing in the ASSEMBLER language results 
from the fact that an element of information on occupies 
one byte, and a minimally addressable field from the problem- 
oriented language UMfl) FORTRAN occupies two bytes. The program 
admits an increment of the functions being executed both at the 
stage of file formation, as well as after the storage procedure 
is completed. 

The results of the operation of the program can be used to 
determine and compare the statistical characteristics of various 
fragments of an image, and also for qualitative and quantitative 
analysis of materials obtained from remote sensing of the earth’s 
surface . 
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2 . Description of the functions Implemented In the program 75 

We shall call a set of elements from several neighboring 
lines which includes elements with the same number of columns an 
oriented image fragment . The program makes it possible to pro- 
cess an arbitrarily chosen oriented fragment of any magnitude, 
up to the entire image. 

In the version of the program to be described the following 
functions are implemented: 

R - executes successively digital read-out of the values 
of the signal in all spectral ranges for all lines of 
the given fragment; 

C - constructs a chart of the line-to-line change in the 
mathematical expectation of the color vectors for the 
elements which are between the extreme columns of the 
given fragment ; 

S - interprets graphically the change in the color vector 
of each element along the lines of the given oriented 
fragment, beginning with its first line. 

M - computes the vector value of the mathematical expecta- 
tion of the signal for the entire given fragment; 

r - executes digital read-out of the files for the frequen- 
cies of the intensities of the fragment in all spectral 
ranges ; 

x - computes the components of the vectors for the mean value 
M* and the dispersion of the signal in the £th 

channel in the given range involved in the signal change; 
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P - graphs the frequencies *<{ and the relative frequencies 
! r <-\ Involved in the intensities of the elements of a 
fragment in any given channel; 

G - prepares a file with elements of the type , /6 

representing a line of a fragment or its part; 

B - computes the self-covariational and self-correlational 
functions of the line. 

The program described consists of three sections of ASSEMBLER 

t 

tiiVTRtA ,*«■»«* .MfcctHtt) [file]) and several modules of the type , 

written in the nofl FORTRAN (Pig. 1). In order for the program 
to begin operating, initiation on the part of the calling program 
is necessary in whose function there is also a storage reserva- 
tion for the operating file . The scale of the file is 

determined by the number of spectral ranges. It must be deter- 
mined for the functions by the formula ^ss«kchm< , and for 

the function G — the formula is the length of the part of the line 
x kch*h t where is the number of channels of video infor- 

mation. The file *** for the functions X,F must be described 
as , and the file for G, as |**re,cM 2 . Prom 

the operator control is transferred to the 

section -tetTHM which reads the M, , searches for and isolates 
the given fragments, implements the functions R, C, S, M, r, and 
also prepares the files used by the algorithms of the functions 
P, B, X and G. The sequence of operations for the blocks and 
modules te»tam (Pig. 2) is determined by the symbol for choosing 
a function which together with the coordinates of the fragment 
being processed must be printed on the control punched card. 

The coordinates of a fragment are the number of the initial and 
final lines c*P'* and ***** and the columns and 

The control parameters of this punched card also include the 
number of spectral ranges (channels) kcha**'' 
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Let us consider in more detail the possibilities represented 
by the program. To Implement all the functions of the program 
after entering the control card, the next line Is read and the 
given Initial number Is compared with the number which Is In the 
first four bytes of the line on . If the number read Is less 
than the given number, then reading proceeds until the numbers 
coincide. If the initially read number is greater than the given 
number, then fictitious reading of a group of lines backward 
occurs (the number of lines is determined by the formula ♦ l)j :» 

and the group is read once In the forward direction in the stan- 
dard mode (checking the numbers). 

When the symbols R,C,S are given, printing on the AttJtf begins 
immediately after the next line of the given oriented fragment 
is read. 

In the case of the function G, the subprogram transforms 
the multi-spectral data of a single line into the form 4NTEdER* 2. 
and returns control to the initiating program for further process- 
ing of the file IflR . In this case the length of the output 
file is determined by the formula (*»«. *!%)«««< (here char- 

acterizes the number of color vectors of the initial line). The 
file created is organized as follows: the first elements 

correspond to the first spectral channel, the second, to the 
following channel, etc. To implement the remaining functions, 
files are accumulated for all the given fragments. 

By means of the function R in one line on an pW/ are printed 
the digital values of the elements of a line in the first spec- 
tral range (but not more than 30 items); in the following line, 
the values of the elements of the next spectral range, etc. If 
the number of elements in a line of the fragment under considera- 
tion exceed 30, then again the * <«* of the lines is printed, 
etc., until all the elements of the line are printed out. The 
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Fig. 1. Hierarchy of the blocks and modules of the program. 






♦Illegible In foreign text. 


Fig. 2. Microenlarged block diagram of ^a »*TES TAM/ 











number of the line and the number of the initial column are 
printed before the beginning of the line corresponding to the 
first spectral range. Then from the nil the following line is 
printed, and the process is repeated until the last line of the 
given fragment has been read. After this, control Is shifted to 
the input block of th e next control card. 

As a consequence of the function S each line of the jtfSiy is /10 
a mapping of a color vector for the given fragment. Since the 
width of the •A.tfiy belt makes it possible to accommodate only 
128 symbols (and the value of an element, as stated above, can 
amount at most to 255 conventional units), the values of all the 
components of a color vector are divisible by 2. Before the 
first multi-channel element of each line of a fragment is printed 
on the Aiflis , the numbers of the mapped line and the initial 
column of the fragment are printed. In the positions determined 
by the values of the components of a color vector, from 1 to 8 
symbols are printed. If the Intensities of the signal in 
various channels coincide, then on the print the symbols are 
superimposed on each other. Thus all the elements of one line 
are printed, then all those of another, etc., until the last 
element of the last line in the given fragment is printed out. 

To facilitate inspection of the material, a coordinate network 
is printed over the graph (through 20 elements along the axis of 
abscissas and 40 conventional units along the axis of ordinates 
of the graph). In the nodes of this network, on the axes of 
coordinates along the perforations are printed the numbers 
designating the number of the color vector in a line, and in the 
perpendicular direction, its intensity. During a single access 
to the function S it is possible to map graphically an arbitrary 
number of lines of a fragment. 

Implementation of the function C after reading through the 
next line in the module ffin* the file: 


8 



( 1 ) 


«t*T« 

i*Sw* 


is stored, and the value 


c (f) . 

*«»• MlfTI 

is printed out on the also: 


( 2 ) 


The range of change of each component c(e) also cannot 
exceed 128 conventional units (i.e., the entire width of the 
jinny belt is used) . Reading the next line of Information calls /II 
the AlSl? to print one line, analogous to the function S. The 
number of a line is printed only on the axis of abscissas, and 
before the entire graph is printed, the number of the element is 
printed with which averaging begins. The number of elements of 
a line which can be averaged is arbitrary, and the number of 
lines is limited only by common sense in the consumption of paper. 


As regards the function M the U23 print out the text "the 
mathematical expectation of a color vector of the fragment", and 
in succession all the of the components of the vector for 

the mathematical expectation §* are printed. These components 
are found by means of the formula 


£*• — — ■ 
9 ****** . i 






mns jtuf 


(3) 


The results of executing the function r Is the digital 
mapping of the file , prepared by the module . This 

file contains the frequencies with which particular values of 
each component of all the color vectors involved in a given 
fragment appear. A single line of such a file contains 255 
numbers, corresponding to the entire intensity range for each 
channel. But during printing on the ajuw the quantity of these 
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numbers may be reduced by preliminary averaging with any multi- 
plicity mat , The value of the multiplicity is determined by 
the positions on the control card. 


r * 

Each element % 
means of the formula 


of the file being mapped is calculated by 


r* 
' < 


« 






( 4 ) 


The prepared file is printed on the A!Uiy , channel by channel, 
and the spectral channels are separated from each other by a line 
consisting entirely of the letters . For the 1th channel, 

the values of the file if, are printed successively in a line by 
the aWF , and for each element *. If the number of file elements /12 
exceeds 30, then the remainder are transferred to the following 
line, etc., until all 256 mVt values of the file are 

withdrawn. To make it easier for the investigator to study the 
distribution of the intensities, the values of the intensities 
corresponding to the first element J of each line are printed in 
the first four positions of that line. They are separated from 
the digital values of the frequencies by the symbol % 


In order to Implement the function X (calculation of the 
mathematical expectation and the dispersion) we begin by storing 
the file w* which determines the frequency with which the 
intensities of the corresponding components for the color vec- 
tors appear. Then in the subprogram kmh, we compute the 

components of the vectors of the first and second moments, using 
the formulas 




t 


£ 


f •!* 

f £ 

t* • 


l 


t-k 


t 


r 

«*• 





(5) 


•Illegible in foreign text. 
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and 


k » t*f¥* «• uf 

'SferC' 1 *-) 


( 6 ) 


Here "E denotes the integral part of the given expression, 
and fc is the number of the channel. 

When calculating the characteristics of the change in the 
signal in a selected part of the intensity range, rather than in 
the entire range, we take as the lower and upper limits of 
summation, the minimum and maximum values of the signals for the 
portion of the intensity under consideration. The subprogram 
vxnm uses the following additional control parameters: 


wcham - the number of the spectral range; 

(n*«v»i,HA*vAL - the minimum and maximum values of the 

intensities of the chosen signal change 
range . 


The AfflU prints the line in which appear the number of the 
spectral range of the channel, the limiting values of the signal, 
the value of the mathematical expectation and the dispersion 
(together with the explanatory text). 

The result of executing the function F is the files of /13 

histograms which make it possible to analyze the probability 
density distribution functions of the signal with respect to the 
intensities . 


For the width of the AiW belt, the program permit 2 the 
output of one, two, three and four histograms. The printing error 
is determined by the width of the field removed for the histo- 
grams. Just as in the case of the function X, it is possible to 
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put out a selected part of the signal Intensity range, rather 
than all of It, and to construct histograms with an average 
relative to a given parameter. Here as the quantum for a his- 
togram we take the average value with respect to several quanta, 
the number of which is given by the averaging parameter r s«c* . 
During operation of the algorithm for the function P, the module 
, after forming the histogram files, transfers control 
to the subprogram »*«***• which is the control program for print- 
ing histograms. This subprogram uses the following control 
parameters which assign the mode for constructing graphs: 

NGIST - the number of histograms for the AW belt width; 

NNNNC - the number of spectral ranges which must be mapped 
on the AW. simultaneously according to NGIST; 

NBEG C , NEND e - the beginning and ending quanta of the mapping 

I e - i , . . kch vt) 


ISRED - the averaging parameter; 


VIDG - the type of histrograms (the frequency of impact 
(A) or the relative frequency (F)). 


When constructing histograms reflecting the frequencies , 

the ij-ro quantum represents the number of elements of the given 
fragment with an intensity of i units. The quantum o£ in a 
histogram of relative frequencies is calculated by normalizing 
the quanta * ; l so that 


o; 



M. 


( 7 ) 


The aioiy prints the number of the channel, the greatest frequency /l4 
in the file and the number of this quantum, and also a digit on 
the axis of ordinates. The intensities are printed along the 
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axis of abscissas (the conventional unit of intensity is a 
quantum of a histogram) . The digital network is laid off on the 
axis of ordinates. The frequency in the chosen quantum can be 
calc lated by multiplying an ordinate digit by the number of pad 
symbols, '* % w . For example, in channel 1 (cf. Russian P. 23) 
on the histogram constructed by means of the parameter VIDG*F 
an Intensity of 88 has the frequency 9 x 0. 00*131 3 ■ 0.03882. 


As the result of executing the function B (which determines 
the self-covariant and the self-correlational functions of a line 
of a fragment or its part) the line with number ^ which has 
been read forms a file with elements (Cf. function 

G) and control shifts to the module VINFG, written in FORTRAN. 

The following are computed: the average signal ** for a given 

portion of the line (cf. formula 3) and the values of the self- 
covariant function, using the formula 


u, » 77* c_ 

N .-.I 


The self-correlational function is calculated from the formula 



where *,o, i. ...#*• 

On the control card the user assigns the limiting value of 
the delay K and the portion of the line along which the char- 
acteristics are calculated. The Aliiur prints: the number of the 

line, the number of the column between which the portion of the 
line is located, and the values of the functions for all delays 
from the 0th to the Kth. Control is returned to the module 

for input of the coordinates of the following fragment 
and the new function. 
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In an appendix we present examples illustrating the mapping /15 
of the results obtained by processing several fragments. In a 
separate appendix we give the complete text of the section 
in the ASSEMBLER language, together with the necessary commen- 
taries and the texts of all subprograms, written in qoh, FORTRAN. 

3 . Instructions for working with the program 

The program TECTA is controlled by means of symbols and 
numbers fed in on punched cards . 

Each control card represents a sequence of whole numbers, 
arranged in succession, beginning with the leftmost position on 
the punched card. 

To implement any of the functions, it is necessary to specify 
on a card the coordinates of a fragment, the number of spectral 
ranges and the symbol of the chosen function. The data are 
printed in the following sequence, each of the first five invol- 
ving 4 positions per punched card: 

NBSTR - the number of the beginning line of the fragment; 

NBSTO - the number of the initial column of the fragment; 

KCHAN - the number of spectral ranges; 

NESTR - the number of the final line of the fragment; 

NESTO - the number of the final column of the fragment; 
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* - the symbol indicating that in the next position a func- 
tion will be given (1 position per card); 



VIDFUN - for the choice of a function (one of R,c,S,M,Q, /16 
X.r.F.B). 

In the case of the functions X and F, additional parameters 
must be fed in. 

The function X (for computing the mean value and estimating 
the dispersion of the signal in the given spectral range) uses 
the following parameters: 

NCHAN - the number of the spectral range (1 position) ; 

MINVAL - the minimum intensity in the chosen signal change 
range (3 positions); 

MAXVAL - the maximum intensity (3 positions). 

All these parameters are printed on one punched card. There 
can be any number of such punched cards (with additional para- 
meters). In order to shift to feeding in the coordinates of a 
new fragment and/or specifying another function, it is necessary 
to insert a punched card with a zero in the first position. 

When executing the function X, if the next card in the first 
column is the digit 9» then the program concludes the operation. 

The control punched cards determining the functions F and r 
(graphical or digital output of histograms) have a format differ- 
ent from that of other punched cards: in position 21 is printed 

the averaging parameter ISRED (for F or KPAT) (for r) (each 
occupying 4 columns), rather than the service symbol Next 

comes the symbol *, characterizing the conclusion of the numbers; 
then the symbol of the chosen function (F or r), and, in the 
case of the function F, the additional symbol of VIDG, which 
determines the form of the histogram. This is the symbol 
(A), if the frequencies *<* , and F — in the case of the rela- 
tive frequencies — are calculated and mapped. Next comes the 
punched card which successively determines the initial «*£&* , 


15 



the final quanta of the mapping, beginning with the first 

channel . Each of the parameters and occupy ^ 

* » 

positions. Then comes the punched card on which is given: the /17 

number of histograms on the width of the tiW belt NGIST, and 
the number of the channels relative to which histograms must be 
constructed. After the parameter NGIST comes a space; all the 
parameters of this punched card occupy 1 position. 

The results of processing the fragments presented in the 
Appendix were obtained by means of the pack of control cards 
pictured in Pig. 3 (p. 18). 

b . Suggested ways to modify the program 

In the version described, the initiating program, written 
in ribff FORTRAN, executes only a function reserving storage 
for preparing files accumulated during the operation of indivi- 
dual modules. Prom our standpoint it is advisable to extend 
its capabilities, committing to it the semi-automatic processing 
of results obtained during implementation of other functions. 

For example, it may be necessary for a user to compare the 
characteristics of individual fragments (coefficients of corre- 
lation, moments, etc.) in order to decipher thematically the 
information obtained. Therefore it is necessary to make pro- 
vision for transmitting an initiating program for calculating 
the characteristics of fragments-. They can be processed by means 
of algorithms written in non FORTRAN in the form of subprograms 
called by an initiating program after the next fragment has been 
processed by the module TESTAM (or one of its auxiliary modules — 
the file, ) . 

An even higher degree of automated processing under the 
supervision of an initiating program can be achieved by means of 
the following modification of this program: the section TESTAM 

is altered so that the file of control parameters (coordinates 
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of an oriented fragment and a type of processing function) are 
also transmitted by the initiating program. This makes it 
possible to form these parameters immediately in the initiating 
program (without feeding them in from punched cards). 

One of the simplest methods for modifying an initiating 
program is to enlarge the function Q. This can be achieved by 
adding an arbitrary processing program. (In the version des- 
cribed, a processing subprogram is nominally present). Such a 
subprogram obtains from the module TESTAM a file of lines of the 
type and then it may implement any processing algorithm 

which is conveniently programmable in FORTRAN. 

Besides enlarging the function G, it is advisable to modify 
it so that a transmitted file of the type will con- 

tain information about the entire fragment, and not about an 
individual line of the fragment. This makes it possible to 
process immediately the entire fragment by means of a problem 
oriented language (for example, to apply a binary processing file 
for comparing the textural characteristics of various fragments! 

The set of implemented functions must also be supplemented 
by a function for constructing two-dimensional histograms: the 

plane sections of multidimensional histograms will make it 
possible to analyze the correlational connections of various 
spectral channels. 
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